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ABSTRACT 

'^f^ ' Solar and atmospheric evidences have been established and can be explained by 

\ neutrino masses. Furthermore, other experiments claim a few unconfirmed neu- 

■ trine anomalies. We critically reanalyze the Ozy2/3, LSND and NuTeV anomalies. 

(N ; 

P_il 1. Introduction 
<\ 

■ Solar and atmospheric neutrino data show that lepton flavour is violated, and give 
CN ■ the two established evidences for physics beyond the SM that we have today (the SM 
^ ■ and its first flaw both appeared in 1968). Atmospheric, solar, reactor and u beam 
^ ■ data can be fully explained by neutrino oscillations withM^ 

g: I Am^gl = (2.7 ±0.4)10-3 eV^ sin^ 2^23 = 1-00 ± 0.04, 

g' Am22 = (7.1 ±0.6)10^5 eV^ tan^ ^12 = 0.45 ± 0.06. 

g ! Although the specific dependence on neutrino path- length and energy predicted by 
! oscillations has not yet been fully tested, present data exclude all alternative inter- 
im ! pretations which have been proposed. 

Qh! It is plausible that the physics behind present discoveries is a 3 x 3 Majorana neu- 

r-| , trino mass matrix. This is in fact what one gets adding to the SM Lagrangian higher 

^ ! dimensional operators (which parameterize the low-energy effects of new physics too 



^ ! heavy to be directly probed): 



+ ^ — — + ••• 

Maybe we are observing the first manifestation of a new scale of nature, ~ 
10^^ GeV, similarly to what happened in 1896, when operators suppressed by the 
electroweak scale v = 174 GeV were first seen as radioactivity by Becquerel. LHC 
will directly explore physics at the electroweak scale in ~2008. 

Accessing the neutrino scale could be not so fast. If the theoretical scheme out- 
lined above is true we have seen 4 of the 9 Majorana parameters contained in rriij. 
Planned oscillation experiments seem capable of measuring all 6 Majorana parameters 
which affect oscillations; neutrino-less double-beta (Oi/2/5) experiments could discover 
violation of total lepton number and measure one more parameter. 

Alternatively, neutrino experiments might discover something else which does not 
fit well in the above scheme. Present neutrino data show a few unconfirmed anomalies: 
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Figure 1: Heidelberg- Moscow Qv2[3 data without pulse- shape- analysis. 

1. A reanalysiJ^ of Heidelberg- Moscow^ Oz/2/3 data claims \mee\ ~ eV. 
claims with small 9 and Am^ ~ eV . 

3. NuTe\^ claims v^/iion couplings ~ 1% away from the SM. 

4. The apparent observation of cosmic rays above the GZK cut-off might be related 
to neutrino masses (z^uheJ^cmb Z with rriy ~ eV) or to problems with energy 
calibration. 

We discuss possible interpretations of the first three hints (within the SM and be- 
yond) and their signals. Since this is not established physics, we unavoidably touch 
controversial issues: I try to present what seems true to me without hiding problems 
and signalling the most controversial points. Hopefully future data will lead to a 
definite conclusion, maybe confirming one or more of these anomalies. 



2. Heidelberg-Moscow 

Neutrino masses distort the end-point spectrum of 13 decays. From studies of 
^H — > ^He e Ve the Mainz experiment sets the 95% CL bound niy < 2.2 eV (Troitsk 
has a similar sensitivity . 

Some nuclei cannot /?-decay. For example 32Ge cannot /3-decay to 33AS that is 
heavier, so it PP decays as ^2^6 — > 34Seee Uei^e {Q = 2038.6 keV). Various experi- 
ments have observed this and other analogously rare SM processes. If neutrinos have 
Majorana masses, the L-violating decay 32Ge — > 34Seee has a non zero amplitude 
proportional to rriee, the ee entry of the Majorana neutrino mass matrix. The exper- 
imental signal is two electrons with total kinetic energy equal to the Q value of the 
decay. A reanalysiJ^J of the Heidelberg- Moscow (HM) dat^ shown in fig. claimed 
a (2.2^3.1)cr evidence for Oz/2/3. To properly understand HM data one needs to know 
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that the 0^2(3 signal is a peak at Q = 2038.6 keV with known width, aE ~ 1.6 keV 
given by the energy resolution, emerging over the (3(3 and other backgrounds, which 
are not well known. The evidence iri^ was claimed assuming 

data = (flat background) + (Oz/2/5 peak) 

and distinguishing the two components trough a spectral analysis, restricted to data 
in a small search window around Q with few aE size. The continuous line in fig. E 
shows the resulting best-fit: we obtain a 2.0a evidence. The evidence depends on the 
choice of the window size: we have chosen the window where data mostly look like 
a peak, maximizing the evidence for peak'^. Employing a large window (which would 
be the right choice, if the background were really fiat) the evidence decreases to 0.8cr. 

Furthermore, the HM spectrum contains a few other apparent peaks, some around 
the energies of faint 7 lines of ^^^Bi (a radioactive impurity present in the apparatus). 
Fitting all data using all the information we have 

data = (Bi peaks) + (cte background) + (Oz/2/5 peak) 

we find a 1.2(7 evidence*^ as illustrated by the dashed line in fig. ^ Other reanalyses 
performed along similar lines and precisely described by their authors find Lla-*"^, 
1.46 J^, less than IcP^ (these authors combine HIvP and IGEX^^i ^^^^^^ Bq^-j^ 
experiments use ^^Ge with a similar energy resolution and background level. IGEX 
has about 5 times less statistic and finds a slight deficit of events around the 0^2(3 
Q value, where HM finds a slight excess). In conclusion, the hint for 0u2(3 is not 
statistically significant. HM data would contain an evidence for 01/2(3 if one could 
show that the background around Q is lower than what fig. ^ seems to indicate. 

We now shift topic and study what oscillation data imply on rriee, assuming 3 
neutrinos with Majorana masses. Therefore we rewrite \mee\ in terms of mixing 
angles 9ij, neutrino masses mi 2,3 and Majorana phases a,P 



iTUf, 



cos^ 9i3{mi cos^ 9i2 + m2e*" sin^ Ou) + mse*^ sin^ 6'i3 



''The authors o0 claim that the restriction in the window search is not a critical arbitrary choice, 
in apparent disagreement witl0. The claim i refers to a fit of an average sample of simulated 
data (generated under some assumption), while ref.l^ finds that the restriction in the window search 
turns out to be a critical arbitrary choice when analyzing the real data. 

''The tentative identification of some of the spurious peaks in the HM spectrum with faint ^^^Bi 
7 lines proposed i had been criticized i , because their intensity appeared incompatible with 
the intense ^-'^''Bi 7 lines, clearly present in HM data. This issue has been clarified and the relative 
intensities are now compatible within 2a. The initial estimate has been corrected by a factor of 6 
(as first noticed in footnote 9 of a revised version o^""^', there is a trivial normalization error in the 
published HM data'^') times another ©(few) factor related to pile-up effects (computed assuming 
that •^^"'Bi is located in the copper part of the detector j®. 
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Figure 2: \mee\ range allowed by all oscillation data. 

Using a global fit of all oscillation datsP^, we plot in fig. |2lthe \mee\ range at 90% 
CL as function of the lightest neutrino mass^. The darker regions in fig. |21 show the 
remaining uncertainty in |mee| due to the Majorana phases that we would achieve 
if the present best-fit values of oscillation parameters in eq. (P) were confirmed with 
infinite precision (we are assuming ^13 = 0). 

Combining the bound on \mee\ with solar data which point to less than 

maximal mixing 612 one can derive an interesting bound on the mass m,^ of almost- 
degenerate neutrinoi^ 

< 1.0 (1.5) heV at 90 (99%) CL. (2) 

The factor h ^ 1 parameterizes the uncertainty in the 0i>2f3 nuclear matrix element 
(see sect. 2.1 of ref.^^). Our bound holds under the untested assumption that neu- 
trinos are Majorana particles, and ca n be evaded adding e.g. Dirac neutrino masses. 
Cosmology gives a tighter hmilP^l^, mj, < 0.23 eV at 95% CL, under the untested 
assumption that a minimal inflationary model describes structure formation and can 
be evaded by e.g. compensating u free-streaming with a primordial tilt in the power 
spectrum. 

In the future, the sensitivity of 0i/2/3 experiments to neutrino masses should im- 
prove more significantly than cosmology and /?-decay. Fig.|21shows that planned 0i>2(3 
experiments, which could reach a sensitivity in \mee\ of few meV, should see a signal 
if the spectrum of neutrinos is 'inverted' (i.e. Arn^^ < 0) or 'quasi degenerate' (i.e. 
min rrij > 0.05 eV). 
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model and number of free parameters 
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mainly incompatible with 


main future test 


ideal fit (no known model) 
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AL = 2 decay ft — > ei'^v'e 


6 


12 


Karmen 


TWIST 


3 + 1 : Am^^g^ji^ = Am^gNo 


9 


6 + 9? 


BuGEY + cosmology? 


MiNlBOONE 


3 V and CPT(no Affi^^^) 


10 


15 


KamLAND 


KamLAND 


3 and CPT (no Am^jj^) 


10 


25 


SK atmospheric 


i)^ LBL? 


normal 3 neutrinos 


5 


25 


LSND 


MiNlBOONE 


2 + 2: Am^^^^ii^ = Am2^^ 


9 


30 


SNO 


SNO 


2 + 2: Am2^^,iig = Am2j^ 


9 


50 


SK atmospheric 


LBL 



Table 1: Interpretations of solar, atmospheric and LSND data, ordered according to the quality of 
their global fit. A Ax^ = n'^ roughly signals an incompatibility at n standard deviations. 



3. LSND 

Both the LSNE@l and KarmejJISI 

experiments study z/^j obtained from /i"*" decay 
at rest, that therefore have a well known energy spectrum up to 52.8 MeV. The 
search for z/^ is performed using the detection reaction u^p — > ne~^, that has 

a large cross section. The detectors try to identify both the e"*" and the n (via the 
2.2 MeV 7 emitted when the neutron is captured by a proton). The neutrinos travel 
for L ^ 30 m in LSND and L ^ 17.5 m in Karmen. 

LSND finds an evidence for i/^ u^, that ranges between 3 to 7cr depending on 
how data are analyzed. This happens because LSND has a poor signal/background 
ratio: choosing the selection cuts as i the LSND sample contains 1000 background 
events and less than 100 signal events, distinguished only on a statistical basis. The 
statistical significance of the LSND signal depends on how cuts are chosen, and it 
is crucial that all sources of background have been correctly computed. The main 
backgrounds are cosmic rays and misidentification. The final LSND result P( z/p — » 
z/e) = (2.6 ± 0.8) 10~^ can be explained by oscillations with Am^ > O.leV' (se c ^ 
for the precise range). 

Karmen is cleaner than LSND but has a few times less statistics and shorter base- 
line. Karmen finds 15 events versus an expected background of 15.8 events, excluding 
the part of the (Am^, 9) range suggested by LSND with larger Am^ > 2 eV', while at 
smaller Am^ the longer path-length makes LSND more sensitive than Karmen. The 
LSND anomaly is being tested as z^^ — > z/g by the MiniBoone experiment, which has 
quite different systematics. 

If confirmed, the LSND anomaly will require a significant revision of the standard 
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Figure 3: The LSND region in the 
(sin^ 20LSND, AmLgND/ eV^) plane at 90% and 
99% CL (shaded area), compared with the 90% 
(dashed Hne) and 99% CL (continuous hne) ex- 
clusion bounds from other oscillation data com- 
puted in a 3-1-1 oscillation scheme. 



Figure 4: The SM prediction for {g^^gji) at 
68,99% CL compared with the NuTeV deter- 
mination at 68,90,99% CL (big ellipses). The 
NuTeV central value moves along the PW line 
using different sets of parton distribution func- 
tions that assume s = s and = d". 



picture. In fact, oscillations between the three SM neutrinos are described by two 
independent squared neutrino mass differences, allowing to explain only two of the 
three atmospheric, solar and LSND neutrino anomalies as oscillations. 

One possible global explanation of the three anomalies is that an extra sterile 
neutrino generates one of them. The sterile neutrino can be used to generate ei- 
ther the LSND or the solar or the atmospheric anomaly, or some combination of 
them. Only the first possibility is now compatible with data . According to 
this possibility, named '3-1-1' oscillations in the jargon, oscillations at the 

■^Thc other possibilities would remain strongly disfavoured even if possible reasons for being more 
cautious would apply. R. Foot, hep-ph/0210393 suggests that theoretical errors might have been 
underestimated in analyses of NC enriched SK data which disfavour atmospheric sterile oscillations. 
Furthermore, bounds on the sterile component involved in solar oscillations would be weakened if 
solar models underestimate the ®B flux. H. Pas et al., hep-ph/0209373 suggest that adjusting all 
the small mixing angles allowed by a 4;/ framework might weaken the bounds. 

R. Foot, hep-ph/0210393 suggests that scanarios with i/^^ — > atmospheric oscillations provide 
a global fit of all neutrino data with x^/dof « 291/276, which is acceptable. This is true, but the 
goodness-of-fit (gof) test based on the value of the total is inefficient when dof ^ 1: it may assign 
an acceptable gof probability to a solution which is already excluded. This issue was discussed in 
the context of analyses of solar data i and can be exemplified by recalling that, according to 
global fits of solar and KamLAND datsP^, the LOW solution has been excluded but its naive gof 
is still acceptable (presently it has x^/dof — 89/91). 
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LSND frequency proceed trough z/^ ^ — > z/g, so that the effective z/gZ/^ mixing 
angle at the LSND frequency is 6'lsnd ~ (^esd^s- I^ue to this 'product rule^^ the 
LSND anomaly somewhat conflicts with bounds on des and 9^s from z/g and dis- 
appearance experiments. The situation is quantitatively summarized in tableland 
fig. El The part of the LSND region still marginally compatible with other oscillation 
experiments has Am-LSND ~ ^^"^ sin^ 26'^^ ~ 0.2 (for more details see-'^^ "'■^'). 

Furthermore, cosmology disfavours such 3+1 interpretation of the LSND anomaly 
for two different reasons: because 3+1 oscillations thermalize a fourth sterile neutrino, 
and because it is heavy. Since both issues are still controversial, we now try to 
summarize what cosmology is really telling. 

Massive neutrinos make galaxies less clustered: recent global fits of cosmological 
data find rrii, < 1 eV (the analysis performed by the WMAP team gives the strongest 
bounc^^, criticized bjff^ who find appropriate a more conservative treatment of 
data at small scale and of bias). The bound on mj, gets slightly relaxed if there 
are = A thermalized neutrinodl^. 3+1 oscillations indeed thermalize = 4 
neutrinos before the big-bang nucleosynthesis (BBN) epoch, and = 4 seems not 
compatible with BBN. Assuming standard cosmology, the dominant bound on 
comes from the BBN prediction for the ^He primordial abundancy, using as input 
the baryon asymmetr y ex tracted from CMB data. A global fit of cosmological data 
gives = 2.6 ± 0.3^, which apparently excludes N,^ = A and even disfavours 
A^",^ = 3. One can avoid these conclusions by enlarging the error on the ^He primordial 
abundancy, since its determination is still controversial. In tabled we tried to make 
a quantitative statement, which can be criticized telling either that it underestimates 
or that it overestimates the impact of cosmology. 

When data started disfavouring interpretations of the LSND anomaly based on 
sterile neutrinos, other more exotic solutions were proposed. None of them can fully 
reconcile of data. 

Since solar and atmospheric oscillations have been established in neutrinos but not 
yet in anti-neutrinos, one can try to explain all anomalies assuming a CPT-violating 
neutrino spectrum. The larger Am^ between t'e.^j.r is used to explain the LSND 
anomaly, and the smaller Am^ could be in the atm ospher ic or in the solar range, at 
the price of sacrificing some solaJ^ or atmospheric! ^ 1 '^'^ anti-neutrino data. These 
CPT-violating spectra are now somewhat disfavoured, as quantitatively summarized 
in table ^ KamLAND will soon precisely probe CPT in the solar sector. Doing 
the same in the atmospheric sector would require a dedicated z/^ long-baseline (LBL) 
experiment or a dedicated atmospheric experiment such as Monolith. 

The LSND anomaly might be produced by a non standard (AL = 2!) decay mode 
fL z/gCzJ^. It affects electroweak precision data and the e spectrum in p, decays. 
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The latter signal is under test at TWIST. This /2 decay interpretation of the LSND 
anomaly is disfavoured by Karmen, since unlike oscillations it does not exploit the 
fact that LSND has a longer path-lengthP^ ; a more precise quantitative statement 
than the one in table^will be possible if LSND will analyze their data in this context. 

In conclusion, the recent experimental progress has disfavoured all proposed in- 
terpretations of the LSND anomaly. Some of them are not yet excluded, and will be 
tested by future experiments, as summarized in tabled 



4. NuTeV 

The NuTeV collaboratior0 reported a ~ 3o" anomaly in the NC/CC ratio of deep- 
inelastic muon-neutrino/nucleon scattering. The effective coupling to left-handed 
quarks is found to be about 1% lower than the best fit SM prediction. 

The NuTeV collaboration sent both a and a z/^ beam on an iron target. Scat- 
tering events were detected by a calorimeter. The muon produced in CC events gives 
a long track, while the hadrons in NC events give a short track. NuTeV statistically 
distinguishes NC from CC events putting a cut on the track length. The ratios of 
neutral-current (NC) to charged-current (CC) deep-inelastic neutrino-nucleon scat- 
tering total cross-sections, Ri, and i?p, are free from the uncertainties on the neutrino 
fluxes and contain the most interesting information. We recall the tree-level SM pre- 
diction for these quantities. Including only first generation quarks, for an isoscalar 
target, and to leading order, R^ and R^ are given by 

aiuM uX) , , „ aivM uX) , 1 , 

= v{ = 9l + Rp = -pTT ^ = 9l + -9l, 

a[uj\l jiX) (y[uj\l jiX) r 

where 

_ (t(z/7V ^ fiX) _3q + q 1 
^ ^ (x{uAf fiX) ~ 3q + q ^ 2 
and q = {u + d)/2 and q denote the fraction of the nucleon momentum carried by 
quarks and antiquarks, respectively. In this approximation the single parameter r 
accounts for the uncertain QCD dynamics and, after including various significant but 
'trivial' corrections'^, the NuTeV data can be presented as a measurement of gj^ and 
g]i, the NC/CC ratio of effective z/^g^ and I'^qn couphngs, predicted by the SM to be 

15 5 
gl = -- sin^ + g sin"^ 6'w, = g ^^^^ ^w- 

The big ellipse in fig. HI is the NuTeV result, obtained summing in quadrature statis- 
tical and systematic uncertainties, gi is found to be ~ 3a below its SM prediction. 

"^Cuts, QED and electroweak corrections, iron contains more neutrons than protons, the charm 
threshold,... In principle, only a careful job is needed to include all these effects correctly. 
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The difference of the effective couphngs — g]i ('Paschos-Wolfenstein ratio 
does not depend on r 

_ R,- rR^ ajuAT vX) - o(vM vX) ^ ^ 1 . ^ 

The value of R^, measured at NuTeV is consistent with previous experiments, such 
as CCFr2^. Unhke CCFR, NuTeV has two separate and beams: this is the 
main improvement because, under the above assumptions, allows to get rid of the 
unprecisely known partonic structure of the nucleon using -Rpw The NuTeV value 
of i?pw is ~ 3(7 below its SM prediction. 

Doing precision physics with iron is a delicate task. NLO QCD corrections have 
not taken into account by the NuTeV collaboration. They cancel out in the ideal -Rpw 
observablJ^ and therefore probably cannot explain the NuTeV anomaly (caution is 
needed because what NuTeV really measures is a few 10% different from the total 
cross sections which appear in the ideal PW observable). 

Parton distributions q{x) are extracted from global fits, usually performed under 
two simplifying assumptions: s(x) = s(x) and uP[x) = dP'{x). These approximation 
could fail, at the level of precision reached by NuTeV. In presence of a momentum 
asymmetry g~ = Jq^ x[q{x) — q{x)\dx the ideal PW observable shifts as 

-Rpw = 2 ~ ^^^^ + (EW corrections) + (1.3 + QCD corrections) (m" — — s~). 

Both — (isospin violation) and s~ (strange momentum asymmetry) might pro- 
duce the NuTeV anomaly compatibly with other data. 

Isospin violating effects of order u^—d^ ~ {rriu—md) / Aqcd could reconcile NuTeV 
with the SM compatibly with all other available data. Some detailed computations 
performed replacing QCD with more tractable phenomeno logical models suggest that, 
due to cancellations, isospin- violating effects are somewhat too sma It is not clear 
if a QCD computation would lead to the same conclusion. 

The theoretical situation concerning s~ is similar. Since a nucleon contains 3 
quarks (rather than three antiquarks) one expects that s and s carry comparable (but 
not equal) fractions of the total nucleon momentum. Non perturbative fluctuations 
like p ^ KA are expected to give s harder than s since K is lighter than A. Indeed 

> could explain the NuTeV anomaly. Some model computations suggest that 
s~ is too small^, but again it is not clear how reliable they are. The s~ issue can 
also be addressed relying on experimental data. The global parton fit in'^ found a 
hint for a s~ of the desired sign and magnitude. However, a s~ with the opposite sign 
was suggested by an analysis of charm production data, obtained and analyzed by 
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the NuTeV coUaboratioiP^. This analysis is unrehable for various reasons hsted in 
a note added tcP^. Charm production data probe the strange content of the nucleon 
in a powerful way, but we are not able of reliably guessing what they imply for s~. 
It seems worth performin g a professional analysis, even if this not an easy task. 

Finally, nuclear effects^ could affect Ri, and Rp, but only at low momentum 
transfer and apparently not in a way that allows to reconcile NuTeV data with the 
In part, nuclear effects are automatically included in the NuTeV analysis, 
based on their own parton distributions obtained fitting only iron data. 

In conclusion, testing and eventually excluding such SM 'systematic effects' which 
can produce the NuTeV anomaly seems to be a difficult job*^. Therefore, it is use- 
ful to speculate about possible new physics interpretations which might have cleaner 
signatures. Unfortunately, no particularly compelling new physics with distinctive 
signatures has been found. The problem is that charged lepton couplings agree with 
the SM and have been measured about 10 times more accurately than neutrino cou- 
plings. Proposals which overcome this problem look exotic, while more plausible 
possibilities work only if one deals with constraints in a 'generous' way or introduces 
and fine-tunes enough free parameters. For example, mixing the Z boson with an 
extra Z' boson modifies NC neutrino couplings, but also NC couplings of charged 
leptons {ii and u are unified in the same SU(2)i doublet). New physics that only 
affects the gauge boson propagators cannot fit the NuTeV anomaly due to the same 
constraints. Neutrino oscillations do not work. A reduction of neutrino couplings due 
to a ~ 1% mixing with sterile singlets does not work, because CC neutrino couplings 
have been too precisely tested by fi decay together with precision data. Combinations 
of the above effects with enough unknown parameters can work^. The new physics 
could either be heavy with sizable couplings (so that future colhders should see it) or 
light with small couplings (e.g. a Z' with few GeV mass and negligible mixing with 
the Z). 
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°Using a different terminology, this means that in the meantime in view of 'additional uncertain- 
ties' NuTeV is not inconsistent with the SM. 

^For example, mixing with a sterile singlet together with a heavy higgs fits all data if one 
also assume s th at some extra unspecified new physics affects Mw and therefore ignores this 
measuremenJS^ . A failure of the SM fit of electroweak data would support the case for new physics, 
but at the moment we do not see any convincing problem. 
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